Magnetic phase diagram of a quasi-one-dimensional quantum spin system 
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We propose an analytical ansatz, using which the ordering temperature of a quasi-one-dimensional 
(quasi-lD) antiferromagnetic (AF) system (weakly coupled quantum spin-1/2 chains) in the presence 
of the external magnetic field is calculated. The field dependence of the critical exponents for 
correlation functions of ID subsystems plays a very important role. It determines the region of 
possible re-entrant phase transition, governed by the field. It is shown how the quantum critical 
point between two phases of the ID subsystem, caused by spin-frustrating next-nearest neighbor 
(NNN) and multi-spin ring-like exchanges, affects the field dependence of the ordering temperature. 
Our results qualitatively agree with the features, observed in experiments on quasi-lD AF systems. 

PACS numbers: 75.10.Pq, 75.10.-b, 75.40.-s 



I. INTRODUCTION 

The progress in preparation of quantum spin sub- 
stances with well denned ID subsystems has motivated 
the interest in studies of them during last years. An- 
other reason for the investigation of properties of quasi- 
lD spin systems is the relatively rare possibility of com- 
parison experimental data with results of exact theo- 
ries for many-body models. According to the Mermin- 
Wagner theorem, 1 totally ID spin systems with isotropic 
spin-spin interactions cannot have a magnetic ordering 
at nonzero temperatures. However, for quasi-lD spin 
systems, which ID subsystems have gapless spectrum of 
low-lying excitations, the magnetic susceptibility, specific 
heat, and muon spin relaxation often manifest peculiar- 
ities, characteristic for phase transitions to magnetically 
ordered states. The critical temperature of the ordering 
of a quasi- ID AF Heisenberg spin system was first cal- 
culated in Ref. |2| in the absence of the external magnetic 
field. However, for quasi-lD spin systems, in which (the 
largest) exchange constants along the distinguished direc- 
tion are relatively small (~ 1 — 20 K), the transition to 
the ordered state can be governed by the external mag- 
netic field. Nowadays in low-temperature experiments 
high values of the magnetic field (about 20 T for station- 
ary fields and about 60 T for pulse fields) can be used. 
Therefore, for many quasi-lD spin systems it is possible 
to investigate experimentally how the external magnetic 
field affects the Neel ordering, i.e. to determine the H-T 
phase diagram. 

This is why, the objection of the present study is to 
construct an analytical theory (convenient for compari- 
son with experiments), which has to show how the ex- 
ternal magnetic field affects the magnetic ordering in a 
system of weakly coupled AF Heisenberg spin chains. To 
calculate the ordering temperature of the quasi-lD sys- 
tem we use the mean field approximation for a weak inter- 
chain interactions. In this approach ID subsystems are 
considered as clusters (of infinite size), for the description 
of which we can use non-perturbative results. 3 As a re- 



sult, we propose a relatively simple analytical ansatz for 
the magnetic field dependence for the Neel temperature 
of a quasi-lD Heisenberg AF system. 



II. MEAN-FIELD APPROXIMATION 

Consider a three-dimensional spin-1/2 system with AF 
interactions between spins, which form a hyper-cubic lat- 
tice. In the quasi-lD situation the Hamiltonians of the 
ID subsystems are Heisenberg Hamiltonians of AF spin 
chains 
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where J > is the AF exchange coupling between near- 
est neighbor spins in the chain, H = g/isB, B is the 
magnetic field, g is the g-factor of magnetic ions, and 
Hb is Bohr's magneton. Denote by J' <C J the weak 
inter-chain coupling between spins belonging to different 
ID subsystems of the quasi-lD system. If the system is 
AF-ordered, we can write the magnetization of the n-th 
site of the system as 
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where e x<z are the unit vectors in the x- or z directions, 
M is the average magnetization, and tun is the staggered 
magnetization in the direction, perpendicular to the ex- 
ternal field (the order parameter in the considered case). 
The inter-chain interaction can be taken into account in 
the mean field approximation. In that approximation in 
the AF phase we write the Hamiltonian of the total sys- 
tem as 
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where h^ = zJ'itin, and z is the coordination num- 
ber. The order parameter m^ (or h^) has to be de- 
termined self-consistently. The self-consistency equation 
reads m^ = Mn(H,Iin,T), where Mn{H,1in,T) is the 
magnetization per site of the ID subsystem in the ef- 
fective field H — MzJ 1 at the temperature T. In other 
words, the susceptibility of the quasi-lD system can be 
written in the mean field approximation as 
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and the ordering takes place at the values of the temper- 
ature and the field, at which the denominator becomes 
zero. Then the transition temperature to the ordered 
state has to be determined from the equation 



1 = zJ'xn , 

Xn = (dM N (H,h N ,T)/dh N ) h , 



(5) 



Notice that xn is exponentially small for the situation 
with gapped low-energy eigenstates of the spin chain. 
It takes place, e.g., for the Heisenberg spin chain for 
H > H s in the ground state (H s = 23 is the critical value 
of the magnetic field, at which the spin chain undergoes 
a quantum phase transition to the spin-saturated phase). 
In that case weak couplings J' cannot yield a magneti- 
cally ordered state of a quasi-lD system. Therefore, in 
the following we consider only the case with gapless low- 
energy eigenstates of the spin chain. 



III. SUSCEPTIBILITY OF THE 
ONE-DIMENSIONAL SUBSYSTEM 

The non-uniform static susceptibility of the ID sub- 
system at low temperatures can be written as 



Xa{q,T) 
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(6) 
where q is the wave vector, a = x,y, z, and (...)t denotes 
the thermal average at the temperature T. Asymptotic 
behavior of correlation functions for an integrable spin 
chain for the gapless case can be obtained in the con- 
formal field theory limit r- and it is possible to write the 
staggered part of the correlation function for the trans- 
verse to the magnetic field components in the ground 
state (related to xn) as 



(^)So(O)) » (-1) 



C 



[n 2 - {vt) 2 ]n/ 2 



(7) 



where v is the Fermi velocity of low-energy excitations, 
rj = 1/2Z 2 is the correlation function exponent, Z is 
the dressed charge of low-lying excitations (low-energy 
eigenstates), and C is a non- universal constant.- Asymp- 
totic behavior can be extended for weak nonzero tem- 
peratures using the conformal mapping (n ± vt) — ¥ 



(v/irT) sinh[7rT(n ± vt)/v]. Then, we can calculate sus- 
ceptibilities for q — it (we use the main approximation), 
Fourier transforming of the conformal mapping of Eq. ([7} 
at low temperatures as 
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where B(x,y) = T(x)T(y)/T(x + y) is the Euler's beta 
function. Finally, we obtain the expression for the Neel 
temperature below which the magnetic ordering takes 
place 
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IV. BETHE ANSATZ APPROACH 



Fermi velocity and the critical exponent r\ can be calcu- 
lated exactly using the Bethe ansatz*2 In the ground state 
phase with gapless low-energy eigenstates at H < H s , we 
can write Z = £(A), v — e(A)/2na(A), where £(x) and 
a{x) and s(A) are determined from the solution of the 
Fredholm integral equations of the second kind 
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The boundaries of integrations are related to the value 
of the magnetic field (0 < H < H s for the phase with 
gapless excitations, while for H > H s the spin chain 
is in the spin-saturated phase with gapped excitations) 
via H = 2itJcj{A)/t;(A). Equations ([TOf can be solved 
analytically only in some limiting cases, and numeri- 
cally in other cases. In the absence of interactions be- 
tween z-components of neighboring spins (so-called XY 
model, the Hamiltonian of which can be exactly mapped 
to the one of the non-interacting fcrmion model using 
the Jordan- Wigner transformation), the dressed charge 
is equal to 1. It is also equal to unity for the isotropic 
Heisenberg chain at H = H s , where A = 0. In the 
absence of the magnetic field we have A = oo, and 
the solution of integral equations can be obtained by 
the Fourier transformation, which yields v — it J/2, and 
Z = l/v2. The numerical solution for intermediate val- 
ues of A shows that the dressed charge as a function of H 



grows from \j\/2 to 1 for < H < H Sl see, e.g., Ref. y, 
i.e. J] decreases from 1 to 1/2 in this domain of field 
values. Similarly, the velocity of low-energy excitations 
decreases with the growth of the field from ttJ/2 to zero 
in the domain < H < H s . Numerical solution for the 
Neel temperature was given, e.g., in Ref. 4. 



V. SIMPLE ANALYTIC ANSATZ 

It is not convenient, however, from the viewpoint of 
application of the results for comparisons with experi- 
mental data to use numerical solutions. This is why, we 
propose the simple ansatz for the magnetic field behavior 
of the the velocity v and correlation function exponent 
rj, valid in the interval < H < H s : 

v = ^V[l - (H/H S )}[1 - (H/H s ) + (2H/nJ)} , 
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The non- universal constant is equal to 0.18 at H = 
and near the saturation it behaves approximately as C ~ 
0.18-sA — 2M, where M is the average spin moment per 
site, cf. Ref. y, which leads to the field dependence C ss 
0.18[2(H S - H)/2nH s } 1 / 4 in the vicinity of the critical 
saturation point H s . Our ansatz is exact at the points 
H = and H = H s and in the vicinity of H = H s . The 
main deviations of our ansatz from exact results take 
place for intermediate field values, between zero and H s . 
In the above expression for the Neel temperature (and 
in the ones for the velocity and the critical exponent) we 
did not take into account logarithmic corrections, which 
exist for the characteristics of the isotropic Heisenberg 
AF spin-1/2 chain near H = 0, see, e.g., Ref. yj. Those 
corrections can be taken into account, which yield for the 
susceptibility 6 



Xn -> Xn 



v /ln(24.27J/T) 
(2tt) 3 / 2 



(13) 



Then for the Neel temperature we can modify our Eq. 
as (cfJi): 



G^r {2t:)- 7 ' a \I^/2{H s - H)/H s \n([48.54irJ/v)} . 

(14) 
Equations (J9j) , (|L2[) , and (TT4|) are the main result of our 
work. 



VI. RESULTS FOR THE MAGNETIC PHASE 
DIAGRAM 

In Fig. [1] we present the Neel temperature of a quasi- 
1D spin-1/2 AF system for z — 4 and J 1 = 0.1J as a 
function of the external magnetic field H , i.e. the H — T 
phase diagram of the system (our results qualitatively 
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FIG. 1: The H — T phase diagram of the quasi-lD Heisenberg 
spin-1/2 system with AF interactions between nearest neigh- 
bors (the value of the spin saturation field H s = 2J). The 
dashed line represents the result without logarithmic correc- 
tions, while the solid one is related to the case with taken into 
account logarithmic corrections. 



agree with the results of numerical calculations for Bethe 
ansatz equations, see, e.g., Ref. |4J). The quasi-lD sys- 
tem is in the magnetically ordered state in the interval of 
fields and temperatures, limited by the line of the second 
order phase transition. Logarithmic corrections do not 
change the qualitative behavior of the Neel temperature 
as a function of the field. However, the values of the 
critical temperatures become smaller due to logarithmic 
corrections. The Neel temperature as a function of the 
external field first grows, and then goes to zero at the 
critical field H = H s . Hence, there exists a (narrow) 
interval of temperatures, at which one can observe a re- 
entrant phase transition. In this domain of temperatures, 
if we enlarge the value of the field, a quasi-lD spin sys- 
tem is first in the paramagnetic short-range phase. Then 
the system undergoes a phase transition to the magnet- 
ically ordered phase, and then, for larger values of the 
field, it returns to the paramagnetic phase. It would 
be interesting to observe such a re-entrant phase tran- 
sition in real quasi- ID spin systems. Such a behavior 
is the consequence of the field dependence of the criti- 
cal exponent. If the critical exponent does not depend 
on the magnetic field (e.g., in the XY chain), the Neel 
temperature as a function of the field only decreases (fol- 
lowing the field dependence of the velocity). Exact cal- 
culation of the Neel temperature for quasi- ID spin sys- 
tem with the Dzyaloshinskii-Moriya (DM) interaction 6 - 
also revealed similar to Fig. 1 behavior, i.e. the maxi- 
mum in the field dependence of the critical temperature. 
Our calculations for spin systems with DM interactions 
or with the "easy-plane" magnetic anisotropy show that 
such magnetically anisotropic interactions produce the 
reduction of the maximum in the field dependence of the 



critical temperature. Notice that the critical tempera- 
ture for weakly coupled XY spin chains is higher than 
for Heisenberg chains due to field-independent exponent. 
Also, if the symmetry of the lattice of the total system is 
lower than hyper-cubic, the ordering can take place not 
at q = 7r, but for some values of q, which depend on the 
lattice structure and relativistic interactions (such a case 
can be analyzed in the random phase approximation). 
In that case the Neel temperature also becomes smaller 
than in the hyper-cubic situation, cf. Ref. |g. 



VII. EFFECT OF NEXT-NEAREST NEIGHBOR 
AND MULTI-SPIN EXCHANGE COUPLINGS 

In real quasi-lD spin systems additional intra-chain 
exchange interactions between NNN spins exist very 
ofteni To take into account such interactions we can 
consider the modified ID Hamiltonian 



H 



NNN 



= J\ y y (s ra -s ra -n)+j2 y y (s n -s w +2)— # 7 „ s^ 



where Ji is the exchange integral for next-nearest neigh- 
bor couplings. For J 2 > such a spin chain reveals a spin 
frustration. Unfortunately, for this Hamiltonian an exact 
solution cannot be obtained analytically for any values 
of Ji.2- Nevertheless, approximate bosonization studies 
and numerical calculations suggest that for J2 > 0.24... J\ 
the spin gap is opened for low-energy excitations^ For 
the above mentioned reasons a system of weakly cou- 
pled chains with gapped excitations cannot be ordered 
magnetically. However, as follows from Ref. 7, despite 
the fact that for most of studied compounds exchange 
constants satisfy the condition J2 > 0.24...Ji, the spin 
gap was not confirmed experimentally. To describe the- 
oretically quasi-lD spin systems with spin frustration 
due to intra-chain interactions without spin gap and 
with a weak inter-chain coupling, we consider another 
model, the Hamiltonian of which is 1-Lnnn with addi- 
tional terms, describing multi-spin ring-like interactions. 
The advantage of that model is its exact integrability: 
The model permits an exact Bethe ansatz solution. Wc 
do not state, naturally, that the model describes all fea- 
tures of the experiments. 7 However, many properties of 
the model are similar to what was observed in Ref. @, 
at least, for this model low-energy eigenstates are gap- 
less. Hence, from this viewpoint, it qualitatively agrees 
with the data of experiments p unlike the model with 
the Hamiltonian T-Lnnn- Multi-spin ring exchange inter- 
actions are often present in oxides of transition metals, 
where a direct exchange between magnetic ions is compli- 
mented by a super-exchange between magnetic ions via 
nonmagnetic ones. 9 The modified Hamiltonian of such 
ID subsystem has the form 

liraod — T~LnNN + Ji / , ((S n — 1 S n +1 ) (S w S ra +2) 
n 

— (S n _iS n +2)(S n S n +i)) . (16) 



Notice that multi-spin interactions are less relevant from 
the renormalization group viewpoint than two-spin inter- 
actions. Quantum properties of the model can be seen 
from the exact solution^ for the parametrization of cou- 
pling constants J\ = J(l — y), J2 = Jy/2, J4 = 2Jy 
for any J and y (in what follows we consider J > 0, 
y > 0). This exactly solvable model, while being for- 
mally less realistic than the model with the Hamilto- 
nian "Hnnn, reveals features, more similar to the proper- 
ties of experimentally studied quasi-lD systems with spin 
frustration. 7 For y = the model describes the Heisen- 
berg spin- 1/2 chain. The ground state of the model de- 
pends on values of the parameter y and an external mag- 
netic field. 10 At T = for large values of the magnetic 
field the model is in the spin-saturated phase, divided 
from other phases by the line of the second order quan- 
tum phase transition. For low values of y and H the 
model is in the phase, which properties are similar to 
the phase of the Heisenberg spin-1/2 chain in a weak 
magnetic field (Luttinger liquid)^ The model is in this 
phase for y < y cr = 4/-7T 2 at H = and for y < y cr (H) for 
nonzero fields. The point y cr is the quantum critical one. 
For y > y cr the model is in an incommensurate phase 
with nonzero spontaneous magnetization at H — 0. Last 
two phases are divided from each other by the line of the 
second order quantum phase transition. Quantum phase 
transitions can be observed in the temperature behav- 
ior of thermodynamic characteristics of the model, like 
the magnetic susceptibility and the specific heat, that 
were also calculated exactly/i - This model permits us to 
know, how NNN interactions (together with the multi- 
spin ones), which cause the quantum phase transition, 
can modify the H — T phase diagram, obtained above 
for the case of only nearest neighbor couplings in chains. 
In this situation the expressions for the velocity and the 
critical exponent can be modified using the substitution 
ttJ —> 7rJ(l — x) (cf. Ref. [Io|), where x = y/y cr - We 
concentrate on the case y < y cr (0 < x < 1) for the ID 
subsystem, which has Luttinger liquid properties. In this 
case a quasi- ID spin system undergoes the transition to 
the AF ordered state. 11 Ordering in the incommensurate 
phase was studied in Ref. [Tl]- The results of our analy- 
sis are presented in Fig. 2. The ordered phase is inside 
the region, limited by the surface, at which the second- 
order phase transition takes place. The maximum in the 
field dependence of the critical temperature, cf. Fig. 1, is 
shifted towards low values of the field with the growth of 
x, i.e. spin-frustrating NNN and multi-spin couplings can 
reduce the domain of temperatures, at which re-entrant 
phase transition can take place. We believe, that while 
the considered model seems less realistic, the mentioned 
feature has the generic nature for quasi-lD spin systems. 

Phase diagrams, similar to the ones, presented in 
Figs. 1 and 2, were obtained experimentally for real 
quasi-lD compounds 12 . Phase H — T diagrams in those 
compounds show maxima of the field dependencies of 
critical temperatures. Namely, the ordering temperature 
in studied quasi-lD copper oxides first increases with the 



those compounds measurements reveal different values of 
the ordering temperatures for different directions of the 
external field, which can be caused by the weak magnetic 
anisotropy of the intra-chain exchange interactions. 



T/J 




FIG. 2: The phase diagram for a quasi-lD spin-1/2 chain 
with spin frustration, caused by nearest, NNN interactions 
and the ring exchange. The Neel temperature is a function of 
the parameter x, which shows how close the quantum critical 
point (caused by spin-frustrating interactions) x — 1 is, and 
the external magnetic field. 



growth of the value of the held, reaches its maximum, and 
then decreases to zero at the value of the field, where the 
spin chain has the spin saturation. Notice, that in one of 



VIII. CONCLUSIONS 

In summary, we have used a simple analytical ansatz 
to calculate the ordering temperature of a quasi-lD sys- 
tem, consisting of weakly interacting quantum spin-1/2 
chains with AF couplings in the presence of the exter- 
nal magnetic field, when the weak inter-chain coupling 
is taken into account in the mean field approximation, 
and the characteristics of spin chains are obtained non- 
perturbatively. Our results show that the field depen- 
dence of the critical exponents for correlation functions 
of ID subsystems plays a very important role. In particu- 
lar, that dependence determines the region of possible re- 
entrant phase transition, governed by the field. We have 
shown also how a quantum critical point between two 
phases of the ID subsystem, caused by spin-frustrating 
NNN and multi-spin ring-like exchanges, affects the field 
dependence of the ordering temperature. Our results 
qualitatively agree with the features, observed in exper- 
iments on quasi-lD AF systems. We expect that our 
results are generic for quasi- ID systems and that they 
can be helpful for experimentalists, who study magnetic 
properties of such systems, especially due to the recent 
progress in obtaining high values of the magnetic field in 
experiments. 
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